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Abstract 

Genomic and genetic analyses have demonstrated that many species contain multiple chemotaxis-like signal transduction 
cascades that likely control processes other than chemotaxis. The Cne^ signal transduction cascade from Rhodospirillum 
centenum is one such example that regulates development of dormant cysts. This Che-like cascade contains two hybrid 
response regulator-histidine kinases, CheAs and CheSs, and a single-domain response regulator CheYs. We demonstrate 
that cheS^ is epistatic to cheA^ and that only CheS3~P can phosphorylate CheY3. We further show that CheA3 derepresses 
cyst formation by phosphorylating a CheSs receiver domain. These results demonstrate that the flow of phosphate as 
defined by the paradigm £ coli chemotaxis cascade does not necessarily hold true for non-chemotactic Che-like signal 
transduction cascades. 

Citation: He K, Marden JN, Quardokus EM, Bauer CE (2013) Phosphate Flow between Hybrid Histidine Kinases CheAs and CheSs Controls Rhodospirillum 
centenum Cyst Formation. PLoS Genet 9(12): el004002. doi:10.1371/journal.pgen.l004002 

Editor: Patrick H. Viollier, University of Geneva Medical School, Switzerland 

Received March 22, 2013; Accepted October 21, 2013; Published December 19, 2013 

Copyright: © 201 3 He et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: Funding for this study was provided by a National Institutes of Health grant GM099703 to CEB. The funders had no role in study design, data collection 
and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 

* E-mail: Bauer@indiana.edu 

n Current address: Cystic Fibrosis/Pulmonary Research and Treatment Center, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States 
of America. 



Introduction 

Rhodospirillum centenum is a photosynthetic member of the 
Azospirillum clade, members of which associate witli root 
rhizospheres in a broad range of plants. These aerobic nitrogen 
fixating organisms are capable of promoting plant growth by the 
donation of both frxed nitrogen and plant hormones [1]. 
Inoculating fields and/or seeds with Azospirillum sp. have signifi- 
cantly enhanced crop yields of a wide diversity of cultivars 
including com and wheat [2,3]. An additional feature of this group 
is the capability of forming metabolicaUy dormant cysts that 
promotes survival during droughts [4]. Encystment involves 
several morphological transitions during which cells round up 
and form a thick outer exopolysaccharide coat termed the exine 
layer [5]. The formation of cysts also correlates with the 
appearance of intracellular poly-P-hydroxybutyrate (PHB) gran- 
ules that are presumably used as energy reserves [6] . Once water 
and nutrients are available, cysts germinate by reforming 
vegetative cells that emerge from the exine coat [5]. 

Azospirillum species are morphologically similar to myxospores 
synthesized by Myxobacteria. Both groups are soil-dwelling, 
Gram-negative proteobacteria that form highly desiccation 
resistant resting cells. In Mjxococcus xanthus a two-component 
system (TCS) comprised of a membrane bound histidine kinase 
(HK) CrdS, which phosphorylates a DNA binding response 
regulator (RR) CrdA to control myxospore development. The 
Che-like Che3 signaling cascade negatively regulates CrdA by 



functioning as a phosphatase [7] . As is the case with Myxobacteria, 
cyst formation in R. centenum also utilizes a novel chemotaxis-Hke 
signal transduction cascade (Chcj) to control the timing of 
development [8]. The R. centenum che-j gene cluster (Figxire 1) is 
comprised of eight genes coding for homologs of CheA (CheAj), 
CheW (CheWg, and CheW3b), CheB (CheBj), CheR (CheRs), a 
methyl-accepting chemorecepter (MCP3) and CheY (CheY3). 
CheAj is a CheA-CheY hybrid (Figure 1) belonging to Class II 
HKs, which include homologs of the E. coli CheA with a conserved 
histidine residue located in a histidine phosphotransfer (Hpt) 
domain rather than a dimerization and hisitidine phosphotransfer 
(DHp) domain found in Class I HKs. In addition to CheA3, the che^j 
cluster also codes for a second HK (CheS^). CheSj has two REC 
domains followed by a PAS (Per, Arnt, Sim) domain and a HWE 
Class I HK domain (Figure 1); however, only one of the CheSj REC 
domains contains a predicted phosphorylatable aspartate (D54 in 
RECl, Figure 1) with the comparable position in the second REC 
being substituted by an alanine (A191 in REC2, Figure 1). 

Clearly the presence of a second HK and two additional 
phosphorylatable REC domains in the R. centenum Chcs cascade 
indicates that the flow of phosphate is more complex in this 
signaling pathway than for the E. coli Che signaling cascade. In the 
classic Escherichia coli chemotaxis model, CheA is tethered to the 
MCP-CheW complex and its autophosphorylation at a conserved 
His in the Hpt domain is enhanced upon repellents binding to MCP 
and inhibited upon binding of attractants. CheA phosphorylates a 
conserved Asp in CheY; phosphorylated CheY in turn binds to the 
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Author Summary 

Bacteria use chemotaxis and chemotaxis-like signal trans- 
duction pathways to quicl<ly sense and adapt to a 
constantly changing environment. The purple photosyn- 
thetic bacterium Rhodospirlllum centenum is able to 
withstand long periods of desiccation by forming meta- 
bolically dormant cyst cells, the development of which is 
regulated by the Ches chemotaxis-like pathway. Using a 
combination of genetic and biochemical approaches, we 
demonstrate that hybrid histidine kinase (HHK) CheAs 
encoded in the che^ gene cluster is essential for cyst 
formation while another HHK CheSs inhibits cyst forma- 
tion. We further show that the appended receiver domains 
of these kinases regulate their respective histidine kinase 
domains and are critical in controlling the timing of cyst 
formation. Finally, we demonstrate that CheAs functions 
upstream of CheS3 and promotes cyst formation by 
phosphorylating CheSs. 



flagellum's rotor causing reversal of flagellar rotation. Similar to the 
smooth-swimming and tumbling phenotypes exhibited in E. coli 
chemotaxis mutants, in-frame deletions of individual che^ genes 
produce distinctly opposing phenotypes [8]. Deletions of cheS-j, 
cheT-j, or c/wBj lead to a hyper-cyst phenotype characterized by 
premature formation of cysts, whereas null mutants oimcpj, cheW^^ 
cheWsi,, cheR-j, or cheA^ produce hypo-cyst strains that are defective 
for cyst development [8] . These genetic studies indicate that CheSj 
and CheY3 may constitute cognate partners in a TCS that 
suppresses encystment, and that CheA^ either inhibits phosphor- 
ylation of the CheS.H-CheYg TCS or is part of a separate pathway. 
Here we report that CheY3 indeed accepts phosphates from CheSj 
and not CheAa, and that CheAj derepresses cyst formation by 
phosphorylating the RECl domain of CheSj. 

Results 

Mutations in cheAj, cheSs and cheYs lead to defects in 
cyst formation 

We previously reported that deletions of hybrid histidine kinase 
(HHK) genes c/iM-j and cheS^j lead to opposing defects in the timing 



of cyst formation [8] . Specifically, a deletion of cheA;j resulted in 
severely defective encystment, whUe a deletion of cheS^j resulted in 
enhanced encystment. We also observed that a cheT^ null mutation 
is indistinguishable from the hypercyst phenotype exhibited by a 
null mutation of cheS^. In order to further probe the importance of 
the linked CheAj and CheSj REC domains we introduced alanine 
substitutions at the predicted Asp sites of phosphorylation and 
recombined these mutations into the native R. centenum chromo- 
somal loci (Figure 1). Mutated strains were subsequently assayed 
for cyst development by growth on either nutrient-rich CENS 
medium that promotes vegetative growth or on cyst-inducing 
CENBA medium. Phase contrast microscopy was then used to 
visually assess cyst production coupled with flow cytometry 
quantitation of vegetative/ cyst ceU populations (Figure 2). 

As observed in previous studies, growth of wild type cells in 
CENS medium visibly leads to >99% vegetative cells (Figure 2A), 
whereas growth in CENBA medium produces large cyst clusters 
(Figwe 2B). Separation of individual vegetative ceUs from cyst 
clusters using flow cytometry indicates that the large population of 
vegetative cells present in CENS medium form a tight pattern near 
the origin of a side scatter (SSC) versus forward scatter (FSC) flow 
cytometry plot (Figure SI). In contrast, wild type ceUs grown in 
CENBA medium, which microscopically have a large number of 
cysts clusters, shows a distinct "comet tail" comprised of larger cyst 
ceUs that separate from the tight clustering of smaller single 
vegetative ceUs during flow cytometry (Figure SI). The tight 
clustering of vegetative cells is indicative of a high degree of 
uniformity of ceU size (~1 |J,m) [9] and internal complexity 
whereas the "comet tail" distribution of the cyst ceU population 
shows that there is a wider distribution of sizes (2-8 |J.m) [10] 
present with varying internal complexity due in part to varying 
numbers and sizes of large PHB storage granules inside cysts 
[10,1 1]. Because each cyst cluster typically contains 2 to 6 cells, the 
number of cyst ceUs is significantly higher (estimated to be ~ 4-fold 
higher) than what is measured by flow cytometry quantitation of 
cyst clusters. 

Flow cytometry analysis of wild type cells grown on cyst 
inducing CENBA medium show that ~10% of the cell culture can 
be separated from the vegetative cell population as larger cyst 
clusters (Figwe 2B). In contrast, growth of the AcheA^ mutant in 
cyst inducing CENBA shows a two-log reduction in cyst formation 
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Figure 1 . Gene arrangement of the R. centenum ches cluster and domain organizations of CKieAs, CheSa, and CheYs. Arrow length Is 
proportional to gene length. Abbreviations: REC, receiver domain; PAS, Per, Arnt, Sim domains; HWE_HK, HWE superfamily of histidine kinases; Hpt, 
histidine phosphotransfer domain; CA, catalytic and ATP-binding domain. Conserved histidine and aspartate residues as putative phosphorylation 
sites are denoted for each protein. The start and end amino acid positions of the receiver domains as well as those of the full proteins are also labeled 
according to the prediction by SMART [48]. 
doi:1 0.1 371/journal.pgen.1 004002.g001 
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Figure 2. Characterization of cheS^ cheA^ and cheY^ mutants. The encystment phenotypes of 11 strains including wild type and single, 
double, and point mutants of cheSj, cheA^ and cheY^ were measured qualitatively by phase contrast microscopy and quantitatively by flow 
cytometry. (A) Growth on nutrient-rich CENS medium reveals hyper-cyst strains that overproduce cysts relative to wild type. (B) Growth on nutrient- 
limiting CENBA medium identifies hypo-cyst strains that under-produce cyst cells relative to wild type. Error bars in the bar graphs represent standard 
deviation obtained from two biological replicates. * (p<0.05), ** (p<0.01) when compared to the wild type (wt) strain in an unpaired f-test. 
doi:10.1371/journal.pgen.1004002.g002 



(Figure 2B) to a level that is comparable with that of wild type cells 
growth in vegetative CENS medium (Figure 2A). Not surprisingly, 
the tfe.4}:H49A HK mutant resembles a AcheA^ mutant, as this 
strain also contains a large predominance of vegetative cells 
irrespective of growth on nutrient-rich CENS or cyst-inducing 
CENBA medium. Interestingly, the cheAf.D663A REC mutant 
exhibits an opposing phenotype in that it forms large numbers of 
cysts in both CENS and CENBA growth media (Figure 2). Indeed 
the level of cyst production by the c&^jiDGGSA REC mutant 
exceeds that of wild type cells grown in CENBA. 

The cyst deficient phenotypes exhibited by the AcheA-^ and 
cheA-j:li4:9A mutants are markedly contrasted by the AcheS^ and 
cheSf.ii4:53A HK mutant strains that produce cysts in both CENS 
and CENBA medium. Interestingly, similar to what was observed 
in the CheAa HK and REC mutant strains, the c&5}:D54A REC 
mutant exhibits a cyst defective phenotype that is opposite of the 
hypercyst phenotype exhibited by the AcheSs and cheS-j:Yi4:b'3A HK 
mutant strains (Figure 2). The opposing encystment phenotypes 
produced by the cheA:^ and cheS^^ HK and REC domain mutations 
indicates that the REC domains have regulatory control over the 
linked HK domains in both kinases. Similar to the AcheS:^ and 
eyie5'5:H453A mutants, both the AcAeT-,, and c/zer,:D64A mutants 



produced cyst cells when grown in both vegetative CENS and cyst 
inducing CENBA growth media (Figure 2). 

Finally, to determine the hierarchy of CheAs and CheSg within 
the Che3 signaling cascade, we constructed AcheAjAcheS^j and 
AcheA-jAcheTs double mutants and assayed for encystment. These 
double mutations resulted in hyper-cyst strains that resemble the 
AcheSy and AcheT-j phenotypes (Figure 2), suggesting that CheAg 
functions upstream of CheS^ and CheYg in this developmental 
signaling pathway. 

Divalent metal cation dependencies to address 
intramolecular phosphate flow within CheAj 

HHKs are generally able to undergo four reactions in the 
presence of ATP and divalent metal cations: (1) autophosphory- 
lation, where the conserved His residue within the HK domain is 
phosphorylated by the adjacent catalytic and ATP-binding 
domain (CA) using ATP as a substrate; (2) autodephosphorylation 
of the phospho-His residue within the HK domain; (3) phospho- 
transfer, where the REC domain dephosphorylates phospho-His 
and transfers the phosphate to its conserved Asp; and (4) 
autodephosphorylation of the phospho-Asp residue within the 
REC domain to yield inorganic phosphates (Pi) (Figure S2). In 
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addition, phosphoryl group transfer from a response regulator 
back to its cognate HK is also possible. This reverse reaction has 
been observed in the EnvZ-OmpR TCS [12] as well as in 
phosphorelay systems involving a Hpt domain where the forward 
phosphorylation reaction (Hisl^Aspl^His2^Asp2) is partially 
reversible (Asp2^His2^Aspl— *Pi) [13]. In the presence of ATP, 
HHKs may therefore exist as a mixture of four different 
phosphorylation states as illustrated in Figure 3A: unphosphory- 
lated, His-phosphorylated (His~P), Asp-phosphorylated (Asp~P), 
and His-and-Asp-phosphorylated (His~P/Asp~P). 

In order to characterize potential phosphorylation states of wild 
type CheA^ and CheS^, we isolated CheAa and CheSg with 
hexahistidine tags at their N-termrni and performed in vitro 
phosphorylation assays. In early experiments we observed little 
radioactive labeling on CheAj with [y- P] ATP in buffers 
containing Na"*" and Mg^"*", which made it difficult to biochemically 
characterize CheAj. Earlier studies showed that potassium but not 
sodium stimulates autophosphorylation of E. coli Che A [14]. 
Additionally, the Salmonella tjphimurium CheY~P autodephosphor- 
ylates at a high rate in the presence of Mg leading to a low 
amount of '^^P protein labeling, whereas in the presence of C??^ 
autodephosphorylation is impeded leading to a high level of ^^P 
labeling [15]. To test whether different metal ions affected HHK 
phosphorylation, we performed kinase assays on wild type CheAj, 
CheSs and on CheA^, CheSj REC domain mutants in 14 buffers 
containing 25 mM Tris pH 7.5 and varying in 100 mM mono- 
valent and 6 mM total divalent salt compositions (Table SI, 
Buffers 1-14). As shown in Figwe 3B, CheAg exhibited nearly 
undetectable labeling in Buffers 1 and 3-7, all of which contain 
NaCl as a monovalent salt. CheAs labeled considerably better in 




Unphosphorylated His~P 



all K^-containing buffers with maximum labeling observed in 
Buffer 9 containing Ca^"*" as the sole divalent ion. When D663 was 
replaced with an alanine, '^'^P labeling was greatly improved in 
nearly all buffer conditions (Figure 3B). The enhanced labeling of 
CheA3:D663A compared with wild type CheAg suggests that the 
N-terminal HK domain transfers the phosphate to D663 in the 
REC domain, which subsequently undergoes rapid autodepho- 
sphorylation. The D663A REC domain mutation would thus 
effectively trap the phosphate at H49, thereby allowing increased 
accumulation of phosphates. Regarding the enhanced phosphor- 
ylation of wild type CheAa observed in Buffer 9, we propose that 
phosphate is captured at both H49 and D663 residues due to Ca^"*" 
mediated inhibition of receiver domain autodephosphorylation. 
This conclusion is further supported by acid-base stability assays 
described below. Unlike CheAg, CheSs shows no particular metal 
ion preference (Figure 3B). CheS3:D54A exhibits much lower '^P 
incorporation in Buffers 2 and 9, which contain Ca^""" as the only 
divalent metal ion. 

CheAs undergoes intramolecular phosphotransfer 
between HK and REC domains 

HHKs are found in most bacterial genomes [16] with the role of 
the linked REC domain not well established in most cases. 
However, in several studies it has been shown that the HK domain 
favors intramolecular phosphotransfer to the linked REC domain 
[17-19]. We tested whether intramolecular phosphotransfer 
occurs in CheAj and CheSs by determining the phosphorylation 
states of the HK and REC domains. To capture the His~P, 
Asp~P, and His~P/Asp~P forms of these phospho-kinases, we 
used an acid-base stability assay based on differential pH sensitivity 
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Figure 3. Metal ion dependent phospKiorylation of CKieAs and CheSs. (A) Four possible phosphorylation states of phosphorylated hybrid 
histldlne kinases (HHKs). (B) Metal cation dependencies of phosphorylation of Isolated HHKs CheSs and CheAs and their receiver domain mutants. 
dol:1 0.1 371/journal.pgen.1 004002.g003 
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of His and Asp phosphorylated residues. Specifically, His~P 
(Figure 3A-2) bonds are labile in acidic conditions but stable in 
basic conditions [20] while acylphosphates like Asp~P (Figure 3A- 
3) are both acid- and base-labile [21]. 

In this experiment, we phosphorylated CheA^ and CheS^ in 
Buffer 9 (containing Ca^""" as the only divalent cation) for 30 min, 
denatured the phospho-proteins with SDS and treated samples 
with Tris buffer, HCl, or NaOH. Samples were then assayed for 
^■'P-labeling by SDS-PAGE, with the assumption being that 
phosphorylation is preserved in a bufifered solution with a 
physiological pH (Tris pH 7.5) and thus would represent 100% 
phosphorylation of the kinases before acid or base treatment. 

In Figure 4A we show that ~50% of wild type CheA3~P was 
hydrolyzed by exposure to 0.1 M HCl and that it increased to 
~90% hydrolysis by exposure to 1 M HCl. This is contrasted by 
>90% hydrolysis of phosphate observed with the mutant 
(CheA3:D663A~P) in both low and high HCl concentrations. 
The different stability profiles of the wild type CheA^ and D663A 
mutant suggest that Asp~P likely exists in the wild type CheA3~P. 
This is confirmed by treatment with NaOH, which dephosphor- 
ylates only Asp~P. In this case nearly 100% CheA3:D663A~P 
withstood high pH while the phosphate on wild type CheA^ is 
extremely labile (Figure 4A). This demonstrates that CheAj:- 
D663A~P is indeed only phosphorylated on a His residue and 
that wild type CheA has the majority (>90%) of its phosphate 
located at D663. Collectively these data suggest that the phosphate 
group flows from the HK domain to the REC domain within wild 
type CheAj. This conclusion is also confirmed by observing direct 
transfer of phosphate from CheA3:D663A~P to a truncated 
version of CheA^ comprised of only the C-terminal receiver 
domain (CheAs-REC) (Figure 4C). Intermolecular phosphoryl 
transfer to CheA3-REC was also detected using the wild type 
CheA3~P as the donor (Figure S3A) that has a linked REC 
domain competing with intermolecular phosphoryl transfer to the 
truncated REC domain. 

Tethered receiver domains in HHKs can either function as an 
intermediate within a multicomponent phosphorylation cascade, 
or as a phosphate sink, removing phosphate from the HK domain 
to impede it from phosphorylating an untethered cognate REC 
domain. We believe the latter is the case with CheA-j as the half- 
Kfe of the phosphate on the CheA3:D663A mutant is nearly 3-fold 
higher (80 min, Table 1) than is observed with wild type CheAs 
(31 min. Table 1). Taken together, it appears that the REC 
domain in CheA3 functions to modulate the phosphorylation state 
of the HK domain by accepting a phosphate that is then rapidly 
lost by hydrolysis. 

In contrast to the acid and base stability of CheA3, CheS3~P is 
only acid-labile (Figure 4B). Furthermore, substitution of the 
predicted D54 phosphorylation site to an alanine in the first REC 
domain does not alter pH sensitivit}'. These results indicate that 
His~P (Figure 3A-2) is the primary autophosphor^iation form of 
CheS3~P. Because CheS3:D54A showed reduced '^'^P incorpora- 
tion in Buffer 9 (Figure 3B), we repeated this assay with CheS3~P 
and CheS3:D54A~P prepared in Buffer 5 (containing both Ca^""" 
and Mg^"") in order to rule out any ion effects imparted upon the 
phosphorylation equilibriums discussed above (Figure S2). We 
observed the same results of high HCl sensitivity and NaOH 
resistance regardless of the buffer conditions (Figure S4). In 
agreement with this conclusion, no phosphoryl transfer was 
detected from CheS3~P to a truncated version of CheS3 
comprised of only the N-terminal CheS3-RECl domain 
(Figure 4D). Interestingly, despite evidence against CheSj intra- 
molecular phosphoryl transfer, the >4 hour stability of CheS3:- 
D54A~P is substantially greater than the 55 min stability 



observed with CheS3~P (Table 1) suggesting that D54 may 
play a role in promoting autodephosphorylation of the HK 
domain. 

CheSa phosphorylates CheY^ 

Based on the CheA-CheY paradigm from E. coli, we tested the 
ability of CheAs to phosphorylate CheY3. In our assays CheA3~P 
and the more stable CheA3:D663A~P mutant did not exhibit any 

detectable ability to transfer a phosphate to CheY^ (Figure 5A, 5B) 
in Buffer 9. Since the E. coli CheY and other response regulators 
exhibit a wide range of binding affinities to divalent metals (K^ of 
0.4-47 mM under pH 6.0-10.0 have been reported [15,22-24]), 
we also assayed CheA3 phosphorylation of CheYs in Buffers 15-21 
with higher (18 mM) total divalent metal concentrations (Table 
SI). This assay condition also failed to obtain phosphoryl transfer 
from CheA3~P or CheA3:D663A~P to CheY3 (Figure S5). 

In contrast to the hypo-cyst phenotype exhibited by nuU 
mutation of cheA^, null mutations in cheS^ and cheT^ both exhibit 
indistinguishable hyper-cyst phenotypes (Figure 2) indicating that 
CheS3 might be the cognate kinase of CheY3. To test whether 
CheSs can phosphorylate CheY3 we phosphorylated CheSs for 
30 min and then added CheY3. Upon addition of CheY3, rapid 
phosphoryl transfer from CheS3~P to CheY3 was observed within 
30 sec (Figure 5C). We also observed that CheS3:D54A is capable 
of phosphorylating CheY3 (Figure 5D) and that the H453A point 
mutation renders Ch(-S3 unable to autophosphorylate (Figure S6). 
Thus, the phosphoryl group appears to transfer directly from 
H453 from CheS3 to CheY3. We also note that CheY3 appears to 
have a fast autodephosphorylation rate similar to chemotaxis 
CheYs [25-27]. 

CheAs phosphorylates the REC1 domain of CheSs 

Since the REC 1 domain of CheS3 is not phosphorylated by the 
tethered HK domain, we questioned whether CheA3 participates 
in the CheSs pathway by phosphorylating the RECl domain of 
CheSs. We initially performed a phosphotransfer assay using 
CheA3~P as the phospho-donor and did not observe CheS3- 
RECl phosphorylation in Buffer 9 (Figure S3B) or Buffer 15 
(Figure S3C). ^Ve reasoned that it may be difficult to observe an in 
vitro intermolecular transfer of phosphate from CheAs to CheS^ as 
the intramolei:ular transfer (rom the HK domain of CheAj to the 
tethered REC domain of CheA3 may outcompete this reaction. 
We therefore repeated the assay using the CheA3:D663A mutant 
as the tethered mutated REC domain would not compete with this 
intermolecular transfer. As shown in Figure 5E, CheA3:D663A 
does indeed transfer a phosphate to the CheS3-RECl domain. 
This transfer from CheA3 to CheS3-RECl also demonstrates a 
level of specificity typically exhibited between cognate HK-RR 
partners, as phosphate does not ffow from CheS3 to CheA3-REC 
(Figure 5F, Figure S3D). 

As shown in Figure 2, a D54A mutation in the CheSj RECl 
domain that would be unable to accept a phosphate in the REC 
domain exhibits a cyst deficient hypo-cyst phenotype. This is 
opposite of the hyper-cyst phenotype exhibited by a H453A 
mutation (Figure 2) that would disrupt CheSs kinase activity. 
These opposing phenotypes suggest that phosphorylation of the 
CheSs REC 1 domain by the HK domain from CheA3 would have 
an inhibitory effect on autophosphorylation of CheSs or a 
stimulating effect on autodephosphorylation of CheS3. This 
conclusion is also supported by genetic and epistasis studies 
which indicates that cheSi; null mutants are hyper-cyst and also 
epistatic to the hypo-cyst phenotype exhibited by cheA-j nuU 
mutants (Figure 2). 
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Figure 4. Identification of intramolecular phosphoryl transfer within CheAs and CheSs. (A) CheA3~P is acid- and all<aline-labile, whereas 
the REC mutant CheA3:D663A~P is acid-labile and base-resistant. (B) Both CheS3~P and its REC mutant CheS3:D54A are acid-labile and alkaline- 
stable. (C) CheA3:D663A~P phosphorylates CheA3-REC truncation protein in Buffer 15 containing and 18 mM Mg^*. (D) Phosphoryl transfer from 
CheS3:D54A~P to CheS3-REC1 truncation protein was not observed in Buffer 15. 
doi:1 0.1 371 /journal.pgen.1 004002.g004 



Discussion 

The R. centenum ches gene cluster encodes a complex 
chemotaxis system of alternative cellular functions 

Chemotaxis and chemotaxis-like signaling pathways represent 
some of the more complex multicomponent signal transduction 
systems present in prokaryotes. A recent bioinformatic analysis of 
450 non-redundant prokaryotic genomes found that 245 con- 
tained at least one chemotaxis-like protein [28]. In these 245 
genomes there are a total of 416 chemotaxis-like systems that 
contain at least an MCP, CheA, and CheW homologs, which 



together are considered a minimum chemotaxis core [28]. 
Together, Che-like signal transduction cascades are known to 
control three classes of function: flagellar motility, type IV pUi- 
based motility (TFP), and alternative cellular functions (ACF) [28]. 
The ACF class comprises approximately 6% of all the identified 
chemotaxis systems, regulating cellular processes such as cell 
development [29,30], biofdm formation [31], exopolysaccharide 
production [32], cell-cell interactions [33,34], and flagellum 
biosynthesis [35]. In fact, most identifiable Che-like signal 
transduction cascades are yet to be genetically disrupted so the 
function of many of these pathways remains to be elucidated. 
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Table 1. Half-lives of phospho-HHKs. 





Proteins 


Half-life (min) 


CheA3~P 


31.2±2.2 


CheA3:D663A~P 


79.6±2.5 


CheSj-P 


55.4±6.1 


CheS3:D54A~P 


>4 hours 



Standard deviations were calculated from two replicate experiments for each 
protein. 

doi:l 0.1 371 /journal.pgen.l 004002.t001 



Chemotaxis systems either exhibit typical chemotaxis architec- 
ture as found in E. coli, or have evolved to include additional 
auxiliary proteins and/ or multi-domain hybrid components. Only 
a few of the more complex Che-like systems containing auxiliary 
proteins have been biochemically and genetically assayed for the 
flow of phosphate among protein components. Consequently, it 



remains unclear whether the CheA-CheY paradigm from E. coli 
will hold true for the many other, and often more complex, Che- 
like cascades from other species. Clearly the results of this study 
indicate that the Chcg cascade from R. centenum differs from this 
paradigm in that CheAs functions to regulate the CheSj-CheYj 
TCS. In some respects this is similar to the Che3 cascade from M. 
xanthus where a CheA homolog controls developmental program 
by acting as a phosphatase to the DNA binding RR CrdA [7]. 

HHKs with appended REC domains are often present in 
organisms that adopt complex life styles such as M. xanthus [36-39] 
and R. centenum [29,40], allowing for added layers of regulation 
within signaling systems. In some cases, intramolecular phosphoryl 
transfer occurs within HHKs. For example, RodK from M. xanthus 
has three REC domains that are all essential for fruiting body 
formation but the HK domain selectively transfers a phosphate to 
its third REC domain [36] . In E. coli, the HK and REC domains 
of RcsC are involved in a HK-^REC^Hpt-^REC phosphorelay, 
which regulates capsular synthesis and swarming [41]. In other 
cases, the receiver domain can either prevent the HK from 
autophosphorylating, presumably by an occluding mechanism 



Time (min) 0 0.5 1 2 5 
CheA~P — I 



CheY~P — 



B 



Time (min) 0 0.5 1 2 5 



CheA -DSeSA-P — 



CheY ~P — 



C Time (min) 0 0.5 1 2 5 



Time (min) 0 0.5 1 2 5 



CheS ~P — 



CheY~P — 



CheS3:D54A~P — 



CheY ~P — 



Time (min) 0 3 10 20 30 40 F 



Time (min) 0 3 10 20 30 40 



CheA -DSeSA-P — 



CheS,-REC1~P — 




CheS •D54A~P 



CheA-REC~P — 




Figure 5. Intermolecular phosphoryl transfer events assayed among CheS^, CheAa, and CheY^. 2-5 |iM CheSs, CheA3, or their REC 
mutant forms were autophosphorylated in 200 |iM ATP for 30 min before 1/10 volumes of 65 mM CheYs or REC domain truncations were added. (A, 
B) Neither CheA,-? nor CheA3:D663A~P are able to phosphorylate CheYs in Buffer 9 containing and 6 mM Ca^+. (C, D) CheSs-P and 
CheS3:D54A~P phosphorylates CheYs within 15 sec of CheYs addition in Buffer 5 containing Na*, 3 mM Ca^* and 3 mM Mg^*. (E) Intermolecular 
phosphoryl transfer analyses of CheA3 and CheS3. (A) CheA3:D663A~P phosphorylates CheS3-REC1 in Buffer 15 containing K* and 18 mM Mg^*. (F) 
CheS3 is unable to phosphorylate CheA3-REC in Buffer 15. 
doi:1 0.1 371/journal.pgen.1 004002.g005 
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[42], or enhance gene expression by interacting witii tiie cognate 
response regulator of the HHK [43]. 

Cysts are a dormant, non-growing state needed for survival in 
poor growth conditions, so the decision to form or impede this 
developmental pathway must involve multiple inputs and check- 
points. In the R. centenum Chca cascade, there are three receivers 
that are capable of accepting a phosphate from two HHKs (CheBs 
is not discussed here since CheB homologs are typically involved in 
MCP modification and not downstream signaling). CheAj and 
CheSs are HHKs containing respective C-terminal and N- 
terminal REC domains whereas CheY^ is a stand-alone receiver 
without an identifiable output domain. The presence of three REC 
domains and two HK domains encoded in this gene cluster 
potentially makes the Ches signaling cascade quite complex with 
the possibility of multiple inputs and check points, which are 
presumably necessary to control the decision to induce cyst 
formation. 

Phosphorylation levels of CheAs and CheSs are 
modulated by their receiver domains, which have direct 
impact on the timing of cyst formation 

We showed that CheA3~P is acid- and base-labile, indicating 
that an intramolecular phosphoryl transfer occurs between the 
tethered HK and REC domains. This transfer is inhibited when 
D663 is substituted with an Ala, giving rise to a His-phosphor- 
ylated CheA3:D663A~P that is stable at high pH. The phosphate 
on CheA3:D663A is much more stable than observed with wild 
type CheAs, indicating that the tethered REC domain likely 
functions as a phosphate sink, attenuating phosphorylation of its 
own HK domain. Fused REC domains serving as phosphate sinks 
are not unprecedented. ClieAY2, a CheA-CheY hybrid in 
Helicobacter pylori has also been shown to use its REC domain as 
a phosphate sink by rapidly dephosphorylating the linked kinase 
domain [27]. 

Unlike CheA.j, the REC 1 domain of CheS j appears to serve a 
different function. ChcS;!~P is acid-labile and base-resistant and 
also does not phosphorylate its receiver truncation (CheSa-REC 1) 
in vitro. This indicates that the CheSs HK domain does not 
phosphorylate its own REC 1 domain. While it is unclear whether 
the CheSj RECl domain directiy interacts with the HK domain, 
it is evident that the RECl domain greatly affects the phosphor- 
ylation state of H453. This is evidenced by the half-life of 
CheS3:D54A~P that is prolonged by many hours relative to wild 
type CheS3~P (Table 1). Furthermore, the CheS3:D54A mutant 
has an opposing in vivo phenotype from a CheS3:H453A mutant 
thereby indicating that the CheS^ RECl domain has regulator^' 
control over phosphorylation of the CheS^ HK domain. Based on 
these results, we propose that D54 stimulates autodephosphoryla- 
tion of the C-terminal HK domain by a mechanism other than 
transferring and accepting phosphates from the CheS3 HK 
domain. Although we do not yet have molecular details on how 
Asp-phosphorylated CheS3 inhibits the HK domain of CheSj, 
genetic and biochemical results clearly suggest that CheA^ 
promotes cyst formation by phosphorylating the RECl domain 
in CheS3. It is likely that Asp-phosphorylated RECl domain 
causes a conformational adoption that either inhibits CheS3 
autophosphorylation or accelerates autodephosphorylation of the 
tethered HK domain. 

The present Chea pathway model involves 

communication between CheAs and CheSa 

The results of this study allow us to establish a working model 
for the Che3 signal transduction cascade in R. centenum (Figure 6). 



Under cyst non-inducing conditions, CheA3 has low basal level of 
kinase activity that directs intramolecular phosphate flow in the 
direction of His^Asp^Pi. The REC 1 domain of CheS3 remains 
unphosphor)'lated so the HK domain of CheS^ operates at a high 
level of activity that effectively transfers phosphoryl groups to 
CheY3. CheY3~P subsequently activates downstream components 
that repress cyst formation (Figure 6A). Upon starvation or 
desiccation (cyst inducing conditions), a signal is sensed by MCP3, 
which fuUy activates the kinase activity of CheAs (Figure 6B). 
Activated CheAs is now able to phosphorylate the RECl domain 
of CheS^j thereby turning off the HK domain of CheSs leading to 
unphosphor)fated CiuY'^ that induces cyst formation (Figure 6B). 

This model also readily explains the opposing phenotypes of the 
c/uAs:D663A and cheS3:D54A REC mutant strains (Figure S7). In 
the cheAfJi663A REC mutant, intramolecular phosphoryl trans- 
fer, which acts as a CheA3 phosphate sink, would be blocked. The 
resulting elevated phosphate concentration at the CheAs HK 
domain would subsequendy lead to elevated phosphoryl transfer 
from CheA3 to the REC 1 domain of CheS3 under both vegetative 
and cyst inducing growth conditions. Constitutive phosphorylation 
of the REC 1 domain of CheS3 by CheA3:D663A would lead to a 
reduction in the HK activity of CheS3 and subsequent reduction in 
phosphorylation of CheYs. Thus, the cheAy.T)66'^A strain should 
have a cyst defective phenotype under all growth conditions, 
which is what is observed (Figure S7A and B). 

For the cheS} RECl mutant, CheA3 is no longer capable of 
phosphorylating the CheS3 RECl domain due to the D54A 
substitution (Figure S7C). Therefore the CheSs HK domain is able 
to autophosphorylate and phosphorylate CheYs under all growth 
conditions. This would result in constitutive repression of cyst 
formation, which is also observed (Figure S7C and D). 

Multiple cyst developmental signal transduction circuits 

require integration 

Even though details of the Che3 phosphorylation cascade have 
been revealed, several features of this pathway still require 
clarification. First, based on the E. coli chemotaxis model, CheA3 
should be activated by an extracellular signal received by MCP3, 
the nature of which is currentiy unknown. Second, it is unclear 
whether CheSs is regulated only from phosphorylation by CheA3 
or if it also directiy senses changes in metaboUsm during 
encystment via a PAS domain. Third, the outputs and the 
downstream components of the Che3 signal transduction cascade 
remain elusive. One possibility is that CheYs passes its phosphate 
onto unidentified downstream components. 

Also not yet reconciled is how the Chcg pathway is integrated 
with the cGMP signaling in R. centenum. This signaling nucleotide is 
synthesized as cells transition from vegetative growth into the cyst 
developmental phase [44]. While this is a newly identified 
signaling pathway, a cGMP responsive CRP-like transcription 
factor has been identified and is required to induce cyst 
development [44]. How these two seemingly independent 
pathways together control the induction and timing of cyst 
formation constitutes a significant challenge in our understanding 
of this Gram-negative developmental pathway. 

Materials and Methods 

Construction of cheSs and cheAj point mutation suicide 
vectors 

cheSs was PCR amplified with 500 bp of flanking DNA as two 
fragments using wild type cells as template for colony PCR with 
primer pairs listed in Table S2. PCR amplified fragments were 
separately cloned and sequenced in pTOPO. Using a Quikchange 
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©His His(P) 



©His His© ^^^-6»^" 

^ ^^ „^, 

©Asp Asf© 
Pi Pi 



©Asp 
CfieYs 



Cyst Formation 



B Cyst inducing CENBA 




H 

CP, 




©Asp Asf© 
^ ClieAa 
PI Pi 



CheSa 



Asp 
CheYs 



Cyst Formation 



Figure 6. Model for regulation of Che^ signal transduction pathway. (A) In the absence of unknown signals, CheA3 Is deactivated; CheS3 
autophosphorylates and transfers phosphates to Its cognate response regulator CheY3; activated CheY3 then Interacts with downstream components 
to repress cyst formation. (B) In the presence of an unknown signal (denoted by a red star), CheA3 autophosphorylatlon Is activated; 
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His-phosphorylated CheAj constantly transfers the phosphates to its C-terminal REC domain, which serves as a phosphate sink. CheAj—P also 
phosphorylates the REC1 domain of CheSa, inhibiting CheSa kinase activity and CheYj remains unphosphorylated. Cyst formation is therefore 
derepressed without activated CheYs. The thickness of the arrows represents the level of phosphate flow. 
doi:l 0.1 371/journal.pgen.1 004002.g006 



(Stratagene) point mutagenesis kit, the D54A mutation was made 
within the 5' cheS:^ fragment liarboring plasmid, whereas a H453A 
mutation was made in the plasmid harboring the 3' cheSj fragment 
using primers described in Table S2. Suicide vector constructs for 
cheSg containing D54A or H453A mutations were then constructed 
by ligating the appropriate 5' and 3' cheS^ fragments directly into 
pZJD29a using external Bamlil and Xhal sites and were internally 
joined by a Bbsl site common to both fragments. After sequence 
confirmation, plasmids were mated from E. coli SI 7-1 {Xpii) into 
an R. centenum AcheS^ strain [8] . Initial recombinants were selected 
for on CENS''™ and second recombinants with chromosomal 
cheSs point mutants were identified by phenotypic (Gm^/Suc^) 
and colony PGR analyses. Suicide vector constructs for 
cheA:i:li49A, cheAy.DGS'AA, and cheT-,:T)64:A were similarly con- 
structed using point mutagenesis primers detailed in Table S2, 
with cheAg internally ligated using a Clal site and cheT^ cloned as 
one fragment. See Table S3 for a complete list otR. cenknum strains 
used in this study. 

Characterization of cellular morphology and % cyst 
formation by flow cytometry 

Two types of media were used to assay for encystment: CENS 
was used for vegetative growth [45], and CENBA for inducing cyst 
formation [46]. Encystment uninduced cells were prepared by 
overnight growth in CENS at 37°C. Encystment induced cells 
were prepared by washing overnight CENS cultures twice in 
CENBA, subculturing 1:40 into CENBA and then incubating at 
37°C for 3 days. 

For microscopic observations, phase-contrast microscopy was 
performed on a Nikon E800 light microscope equipped with a 
100 X Plan Apo oil objective. For flow cytometry, CENS and 
CENBA cultures were diluted in 40 mM phosphate buffer and 
sonicated briefly (~ 1 sec) at lower power to disaggregate cyst cells. 
All samples were stained in 2 |J.M Syto-9 (Life Technologies/ 
Molecular Probes, Grand Island, NY) for 1.5 hours. Syto-9 is a 
permeant DNA stain that was shown microscopically to penetrate 
both vegetative and cyst cells similarly (data not shown). Initially 
fluorescent calibration beads of 880 nanometers and 10 microns 
were used to set the limits for background. After staining, cells 
were diluted ~1:10-1:20 in 40 mM phosphate buffer just prior to 
running to achieve ~1000 events per second on a Becton 
Dickenson FACS Calibur flow cytometer running CellQuest Pro 
data collection software using an argon laser (488 nm). 100,000 
events were collected per sample with two biological repKcates 
analyzed for each bacterial strain grown in each media. Forward 
and side scatter (SSC vs FSC) were plotted in logarithmic scales. 
Hypo-cyst AcheA^^ and hyper-cyst AcheSy strains were used to 
determine the appropriate gating to use for vegetative cells versus 
cyst cells. Flowjo version 10 (Tree Star, Inc.) was used to analyze 
the data and plot the data for publication. Statistical analysis was 
performed using Prism version 5.0 (GraphPad Software, Inc.). 

Protein overexpression and purification 

Coding regions of CheSs, CheAs, CheYs, and the receiver 
domains of CheAj and CheSj (CheAj-REC and CheSj-REGl) 
were PGR amplified from R. centenum genomic DNA with primers 
listed in Table SI. Gel-purified PGR products were cloned into 
pBluescript SK"*" or pGEM-T, sequenced, then subcloned into the 



Ndel and Xhol sites in vector pET28a. pET28a plasmids for 
overexpression of CheSj and CheA^ point mutants were generated 
using the appropriate pZJD29a vector as template for PGR using 
primers detailed in Table SI. All pET28a constructs were 
transformed into E. coli BL21 Rosetta 2 (DE3) cells (Novagen). 
See Table S3 for a complete list of £. coli strains used in this study. 
For overexpression, overnight cultures of E. coli Rosetta 2 (DE3) 
cells were subcultured 1:100 into 1 L LB medium and shaken at 
37°C to an ODgoo of 0.5. Protein overexpression was induced at 
an isopropyl P-D-l-thiogalactopyranoside concentration of 
0.4 mM and cultures were incubated overnight at 16°C with 
gentle agitation. Cells were pelleted by centrifugation and stored at 
— 80°C until further use. For purification of all proteins, cell pellets 
were resuspended and lysed by ultrasonication in lysis buffer 
(20 mM Tris-HCl pH 7.5, 500 mM NaCl, 25 mM imidazole and 
\0% glycerol). Purification was performed on 1 mL HisTrap HP 
(GE Healthcare) columns using an FPLG system. His-tagged 
proteins were eluted in 20 mM Tris-HCl (pH 7.5) buffers with a 
gradient of 25-500 mM imidazole. Fractions containing purified 
proteins were dialyzed into a storage buffer (25 mM Tris-HGl 
pH 7.5, 100 mM NaCl in 30% glycerol) and stored at -20°C 
until further use. 

Autophosphorylation, phosphotransfer and stability 

assays 

Twenty-one Tris buffers containing common mono- and 
divalent metal ions were used in this study for HHK phosphor- 
ylation and phosphotransfer assays (see the fuU list of buffer 
compositions in Table SI). All kinase reactions and phospho- 
transfers were performed in 0.2 mM final ATP concentration 
except for the half-life determination experiments. AH reactions 
were stopped by addition of 6x SDS-PAGE sample loading 
buffer. All phospho-proteins were separated by SDS-PAGE and 
gels were examined by autoradiography on a Typhoon 9100 
scanner (GE Healthcare) located in the Indiana University 
Physical Biochemistry Instrumentation Facility. 

In the metal ion dependency assays (shown in Figure 3B and 
Figure S3), isolated kinases were diluted in the indicated buffers to 
2-5 |xM. Kinase reactions were initiated by adding 1/20 volume 
of ATP/[y-''''P] ATP mix in 25 mM Tris pH 7.5 and allowed to 
proceed for 30 min at room temperature. For phosphoryl transfer 
to CheYs shown in Figure S3, 1/10 volume of 65 mM CheY3 was 
also added to each reaction mixture at the end of 30 min 
autophosphorylation for another 30 min incubation at room 
temperature. 

In assays assessing intermolecular phosphoryl transfer shown in 
Figure 4C, Figure 5, and Figure 6, ATP mixes and protein 
dilutions were made in same buffer as indicated in the text. 2- 
5 H-M kinases were first phosphorylated in 0.2 mM ATP for 
30 min followed by addition of 1/10 volume of 65 |a,M CheYs or 
receiver domain truncations (CheS^-REC 1 or CheA^-REC). The 
time of receiver addition was set to time 0. Phosphoryl transfer was 
then assessed at various time intervals. 

To determine the half-lives of phosphorylated kinases, 2-5 ^M 
CheAs and CheA3:D663A were pre-autophosphorylated in the 
presence of 10 |iM ATP mix in Buffer 9 for 50 min before passing 
through Bio-Rad Micro Spin 6 chromatography columns to 
remove excess ATP. Dephosphorylation was monitored at room 
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temperature by removing 10 |J,L of the filtrates at various time 
intervals. Phosphorylation of the kinases was quantified using 
ImageJ software by integrating the grayscale density of the 
radioactive bands. % Kinase phosphorylation was plotted over 
300 min and data points were fitted to one phase exponential 
decay using Prism. Half-lives of CheS^ and CheS3:D54A were 
measured with the same protocol with the exception that Buffer 5 
was used in place of Buffer 9. 

Acid-base stability assays 

The phosphorylation state of all CheSa and CheAs variants was 
determined by assaying phosphoprotein stability under acidic or 
basic conditions using a non-filter based assay [47]. Kinases were 
allowed to autophosphorylate at room temperature for 30 min, 
after which phosphoproteins were denatured by adding 0.1 
volume of 20% SDS. Aliquots were withdrawn and mixed with 
equal volumes of 0.1 or 1.0 M Tris pH 7.5, HCl or NaOH and 
incubated for 30 min at 37°C before being neutralized with 1.0 M 
Tris-HCl pH 7.5. Samples were then mixed with 6 x loading dye 
and resolved by SDS-PAGE and assayed for phosphorylation by 
autoradiography. Phosphor\iation of the kinases was quantified 
using ImageJ software by integrating the grayscale density of the 
radioactive bands. 

Supporting Information 

Figure SI Characterization of cheSs, cheA3 and cheY3 mutants 
by flow cytometry. The nutrient-rich CENS medium was used to 
identify hyper-cyst strains and nutrient-limiting CENBA medium 
to identify hypo-cyst strains. SSC, side scatter; FSC, forward 
scatter. 
(TIF) 

Figure S2 Five potential phosphorylation events within (HHKs) 

in the presence of ATP. 

(TIF) 

Figure S3 Phosphor)'! transfer from wild type CheA^ and CheS^ 
to the receiver truncation proteins. (A) CheA3~P phosphorylates 
CheAj-REC in Buffer 9 containing K'^ and 6 mM Ca^+. (B) 
CheA3~P does not show phosphoryl transfer to CheSs-REC 1 in 
Buffer 9. (C) CheA3~P does not show phosphoryl transfer to 
CheS3-RECl in Buffer 15 containing K"^ and 18 mM Mg^"*. (D) 
CheS3~P does not phosphorylate CheA3-REC in Buffer 15. 
(TIF) 

Figure S4 Acid-base stability test on CheS3~P and CheSs: 
D54A~P with phosphorylation performed in Buffer 5 containing 
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conditions, resulting in inactivated CheY3 and therefore dere- 
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CheSs-CheY^ TCS is no longer controlled by CheA3, resulting in 
constitutive repression of cyst formation. 
(TIF) 
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